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ABSTRACT 

Let u(x,t) be the possibly discontinuous entropy solution of a nonlinear scalar conser- 
vation law with smooth initial data. Suppose u e (x, t ) is the solution of an approximate 
viscosity regularization, where e > 0 is the small viscosity amplitude. We show that by 
post-processing the small viscosity approximation u e , we can recover pointwise values of u 
and its derivatives with an error as close to e as desired. 

The analysis relies on the adjoint problem of the forward error equation, which in this 
case amounts to a backward linear transport equation with discon tinuous coefficients. The 
novelty of our approach is to use a (generalized) E-condition of the forward problem in order 
to deduce a W 1,0 ° energy estimate for the discontinuous backward transport equation; this, 
in turn, leads us to e-uniform estimate on moments of the error u e — u. 

Our approach does not ‘follow the characteristics’ and, therefore, applies mutatis mutan- 
dis to other approximate solutions such as E-difference schemes. 
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1. INTRODUCTION 


Consider the scalar, genuinely nonlinear conservation law 

(1.1) ^ t Wx,t)] + -^[f(u(x,t))]= 0. 

It is well-known, e.g. [8], that (1.1) may admit many possible weak solutions. To guarantee 
uniqueness (and in fact .//^-stability), we therefore have to restrict our attention to a subclass 
of possible weak solutions. Namely, we select weak solutions of (1.1) which are realizable as 
small viscosity solutions of 

(1.2) ^[u,(x,*)] + ^[/(u,(x,f))] = e~[Q(U'{x,t))], eQ' j 0. 

With this in mind, we recall that u e (-,t) converges strongly in L l to u(-,t), where u(-,t) 
is the unique, so-called entropy solution of (1.1). The L l convergence rate in this case is 
upper bounded by 

(1.3) \\u e (-,t) -u(-,t)\\ L i < Const. >/e. 

Consult [6], [12] for the discrete analogue of monotone difference schemes, and [14] for spectral 
viscosity approximations. 

We are not satisfied with this error estimate for two related reasons. 

1. Ideally, one would like to recover the entropy solution, u(x,f), within 0(e) error. Al- 
though the estimate (1.3) is sharp, it fails to do so because of the following reason. 

2. The error estimate (1.3) is a global one, while in many practical cases one is interested 
in the local behavior of u(x,t). Consequently, when the error is measured by the 
L 1 -norm, there is a loss of information due to the ‘poor’ resolution of shock (and, in 
the general case, also contact) waves in u(x,t), by the corresponding viscous layers in 

u t (x,t). 

In this paper we are concerned with the local convergence rate of the small viscosity 
solutions, u e (x,t), towards the entropy solution, u(x,t). Assume that initially, at t = 0, the 
initial conditions of the small viscosity problem, u e (-,0), are consistent with smooth initial 
conditions of the conservation law (1-1). Then for any t > 0 we show, in Section 4, that by 
post-processing the small viscosity solution u e (-,t), we can recover pointwise values of the 
(possibly discontinuous) entropy solution u(-,i) and its derivatives, with error as close to e 
as desired. It should be emphasized that our approach does not “follow the characteristics” 
and therefore could be extended to certain discrete approximations for nonlinear conservation 
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laws such as E-difference schemes [10]. Indeed, the present study was originally motivated 
by recent numerical experiments reported in [17]. By post-processing of Spectral Viscosity 
(SV) approximations [16], we were able to recover the pointwise values of discontinuous 
conservative solutions with spectral accuracy, in agreement with the formally spectrally 
small viscosity regularization of the SV method. 

The paper is organized as follows. In Section 2 we study linear transport equations with 
possibly discontinuous coefficients. We show that these discontinuous transport equations 
are well-posed in W l, °° , provided their coefficients are one-sided Lipschitz continuous. 

Such (backward) transport equations arise as the adjoint problems for the forward error 
equation governing the difference u e — u. The coefficients of these backward transport equa- 
tions are indeed upper-sided Lipschitz continuous, in view of Oleinik’s E-condition which 
characterizes the entropy solution. This enables us, in Section 3, to derive e-uniform esti- 
mate on the IV 1 ' 00 -moments of the error u e — u, which in turn, leads to the local recovery of 
the entropy solution discussed in Section 4. Finally, we note that our W 1,0 ° upper-bound on 
the moments of the error provides an independent one-dimensional proof of the usual 0( y/F) 
- Z/ 1 -convergence rate mentioned earlier in (1.3). 

2. LINEAR EQUATIONS WITH DISCONTINUOUS COEFFICIENTS 

In this section we study the linear transport equation 

( 2 . 1 ) + °( X) < )^[^( x > 0 ] = °> 

subject to prescribed initial conditions 

(2.2) <f>(x,t = 0) = <^o(x). 

Assume that the initial data, <f>o(x), are Lipschitz continuous. The standard theory tells us 
that if a(-,f) is sufficiently smooth, say C 1 , then there exists a unique generalized solution, 
of (2.1), (2.2), which remains Lipschitz for all time, 

(2.3) ||<^(-, i)j| L i p < Const T • UMOIIup, -T < t < T. 

The key issue that we address in this section is, roughly speaking, the following question: 
What is the minimal degree of smoothness required from a(-,£) in order to retain the 
Lipschitz continuity of the solution </>(-,£) in (2.1)? 

The answer provided in our next theorem is in the heart of the matter. 

THEOREM 2.1. Consider the linear transport equation (2.1) with Lipschitz continuous 
initial data (2.2). We assume that 
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I. (Uniform Boundedness). o(x,t) is uniformly bounded. 

II. (OSLC). a(x,t) satisfies the following One-Sided Lipschitz Condition. 

(2.4) £-[«(.,«)] 5 M.bf mtl'lO.T], 

Then for t > 0, there exists a unique Lipschitz continuous solution, <f>(x,t), of (2.1), (2.2), 
such that the following estimate holds 

(2.5) Htf-.OIIin.- ^ IIMOIk 1 - ' eM ^\ M(t) = m{r)dT , t > 0. 

Before we turn to the proof of this theorem, a couple of remarks is in order. 

1. Theorem 2.1 allows the coefficient a(-,t ) to be discontinuous ; in fact, increasing jumps 
are permitted. The OSLC assumption (2.4) requires only the decreasing part of a(-,t) 
to be Lipschitz. 

2. Theorem 2.1 asserts that the transport equation (2.1) is well-posed only for positive 
time. In general, linear hyperbolic equation with smooth coefficients is a prototype of 
reversible process. However, since in our case the augmenting OSLC assumption (2.4) 
is irreversible , so is the final conclusion of the theorem. Indeed, simple counterexamples 
can be constructed (see below), which demonstrate that the backward solution <f>{-,t) 
may cease to be Lipschitz in a finite negative time t < 0. 

PROOF. The proof consists of the usual three steps of regularization, a priori energy esti- 
mate and compactness arguments. 

Step 1 . (Regularization). Let ( s (x) = |C(f) be a standard positive C£°-mollifier with unit 
mass. We regularize a(-,f) by spatial convolution, <!$(•, t) = a(-,t) * C&> an< ^ consider the 
regularized equation 

(2.6) J^IMM)] + o. s (x,t)-^-[<f> 6 (x,t)] = 0, *>0 

(2.7) 4>s(x,t = 0) = <j>os{x), <f>os = (f>o* Cs- 

Since we now have a smooth and uniformly bounded coefficient, i), there exists a classical 
C°° solution, <j> s (-,t), of (2.6), (2.7). Since this solution propagates with finite speed, a further 
approximation by truncation (which is omitted) can be used, so that we may restrict our 
attention to the compactly supported case where t)eC^ . Clearly, we have 

(2.8) ||M*»0I k°° ^ ||^o(')IU“- 
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Step 2 . (W 1,0 °- energy estimate). We want to show that ^ is uniformly bounded w.r.t. x,t, 
and 8. To this end we shall carefully iterate on the T p -norms of i). Differentiation of 

(2.6) implies 

or, equivalently 

d d d 

(2.9) + a 6(x,i)-Q^[i>s{*,-t)} = — 7 ^[a s {x,t)]ij) S (x } t). 

Integrating (2.9) against ,p even, over the compact support, we find 

+ J x ^s(x,t)^-[ip p s (x,t)]dx = p -^-[a { (x,t)}rf(x,t)dx, 
and after further integration by parts on the left we arrive at 

(2.10) = (p- ^ J x ~^[as{x,t)]^(x,t)dx. 

We now invoke the OSLO assumption on the coefficient a(-,t). Since (g was chosen as a 
positive mollifier, ag = a * (g satisfies the same OSLO (2.4), namely 

(2.11) ~[ ai (x,t)]<m(t). 

Inserting this into the RHS of (2.10) yields for p even 

(2-12) < (P~ IWOllW-.OIIt- 

Therefore, 

(2.13) \\M;*)\\i* - II - 0)||z,p • eV ww , M(t) = f m^dr, t > 0, 

Jo 

and by letting p f oo we conclude 

( 214 ) II^M-,<)IU» ^ II^MOIU- • e M(i) < ||^o(-)llLi P * e M ^\ t> 0. 

Since a (and likewise ag) is assumed to be uniformly bounded, we can use the equation (2.6) 
to upper-bound the temporal derivative as well, 

(2.15) ||^«-,t)|| i -<||a(.,i)|b-.||*()llL i p-e" W , < > 0. . 

Step 3. (Compactness). By (2.14), (2.15), the uniformly bounded family {<£$} is equicon- 
tinuous - in fact equiLipschitz, and shares a common compact support. Therefore, we can 
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extract a subsequence, still denoted by <f>s{x,t), which converges uniformly to a Lipschitz 
limit function, <f>(x,t). We observe that 

- 1 weak — * in L°°, 

“ fc ^^(-.0 weak- * in Z.°°, 
a*(-,i) -* a(-,<) strongly in L 1 . 

Passing to the limit £ j. 0 in (2.6) we conclude that <^(x, t) is a Lipschitz (generalized) solution 
of the transport equation (2.1), (2.2). Similarly (2.5) follows from (2.8), (2.14) and we are 
done. 

REMARKS. 

1. The H^ 1 * 00 - energy estimate revisited . Consider the 1/ a priori estimates in (2.13). 
The case p — 1 (odd p’s can be justified by further approximation argument which we omit), 
leads to usual BV-contraction of the solution operator associated with (2.6). The proof of 
Theorem 2.1 hinges on the observation that with the help of the OSLO assumption 2.4, one 
can iterate on higher L p norms at the expense of an additional bounded (w.r.t. £) exponent. 
An alternative proof of this essential L°° bound, (2.14), is provided by the following duality 
argument. If A$(x,f) solves the adjoint equation of (2.6) 

A[X«(*.0] + = 0, t < T, 

then A 5 (x,f) = / x A$(£,f)d£ satisfies the backward transport equation 

(2.16) ^[Ae(x,t)\ + a s (x,t)-^[A s (x,t)] = 0, t < T. 

Integrating this against sgnA$(x,t) we obtain — in view of the OSLO assumption (2.4), 

(2.17) ||As(-,t)|| t i > ||A t (-,0)||i. ■ t>0, 

which is the dual estimate of (2.14). 

2. Characteristics . The transport equation (2.1) is governed by the evolution of the 
characteristics 

(2.18) x = a(x,t). 

If a(-,i) is at least Lipschitz continuous, then the ODE (2.18) leads to a ‘nice’ unique 
reversible flow. For our purpose, however, we need less. Namely, in order for to 
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remain Lipschitz for t > 0 . we have to guarantee that while tracing characteristics backward 
in time, different characteristics are not pulled apart. The OSLO assumption is sufficient to 
guarantee that (2.18) generates such a ‘nice’ Lipschitz flow backward in time. The following 
examples demonstrate our point. 

Example 1 . We set 

(2.19) a(x,f) = sgnx. 

Then the Lipschitz continuous solution of (2.1) is given by 

(2.20) <f>(x, t) = <f> 0 ((x - t) + + (x -f t)~). 

Example 2 . We reverse the sign in (2.19), 

(2.20) a(x,t ) = — sgnx. 

Now a(-) has a decreasing jump, and this case is not covered by Theorem 2.1. Indeed, 
different forward characteristic solutions of 

(2.21) x = —sgnx 

may impinge one on each other at finite time. Hence, the backward flow generated by (2.21) is 
ill-posed, and the hyperbolic solution transported along these characteristics cannot remain 
Lipschitz for t > 0. 

Example 3 . We consider the equation 

(2.22) x m (T - x)“, 0 < a < 1, 

which is ‘solved’ backwards (of course, the solution is not unique in this case), starting with 
the final time T > 0. Here, a(x) = (T — x)“ fails to satisfy the OSLC assumption (2.4), 
although the analogue one-sided Holder condition holds. Consequently, the corresponding 
forward transport solution will be only Holder continuous. In contrast to (2.22), let us 
consider the case 

(2.23) x = a(x), a(x) = — (T — x)“, 0 < a < 1. 

Although a(x) is only Holder (but not Lipschitz) continuous, the OSLC assumption (2.4) is 
now fulfilled thanks to the judicious minus sign in (2.23), and the result of Theorem 2.1 is 
valid. 

In the next section we study the error equation associated with the conservation law (1.1) 
and its viscous regularization in (1.2). The dual problem of such an error equation leads to 
a backward transport equation like (2.1). For future reference we therefore state 
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THEOREM 2.2. (The backward transport equation) . Consider the linear transport equa- 


tion 

(2.24) + °( x > x » *)] = °> * ^ T > 

with Lipschitz continuous data prescribed at t = T, 

(2.25) <j>{x,t = T) = <f>(x). 

We assume that 

I. (Uniform Boundedness). a(x,t ) is uniformly bounded. 

II. (OSLC). a(x,t) satisfies the following One-Sided Lipschitz Condition. 

(2.26) L + [o(-,t)] = ess sup M x >*) ~ < m(t), meL l [t,T). 

*#V \ x y ) 

Then fort < T, there exists a unique Lipschitz continuous solution, <p(x, t), of (2.24), (2.25), 
such that the following estimate holds 

(2.27) |W-,*)lk l >~ < IM-)lk‘.~ ■ e M(t) , M{t) = J t m(r)dr, t<T. 


3. A PRIORI ESTIMATE ON MOMENTS OF THE ERROR 

We return to the scalar, genuinely nonlinear (say, strictly convex) conservation law, 

(3.1) 


A[„( x ,<)] + |-[/(u(x,i))]=0, i>0, /"><*> 0. 


It is well-known, [8], [15], that the unique entropy solution of (3.1) is characterized by 
Oleinik’s E-condition [11] 


(3.2) 


u(x,t) -u(j/,t) < 1_ t>Q 

x — y ~ at’ 


Later on we shall need a slightly stronger version of this E-condition — interesting for its own 
sake, which is the content of our next result. 


|[u.(x,t)| + J-[/(«,(*,t))] = c^jWKOM))], C<?' > o. 


THEOREM 3.1. Consider the (possibly degenerate) parabolic regularization of (3.1), 
(3.3) 

Then the following estimate holds 

0 . 4 ) ^™W) >0 - 
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PROOF. Let us first assume that Q' is strictly positive, so that the uniformly parabolic 
equation (3.3) admits a smooth solution. Multiplying (3.3) by (J)'(u e (x, f)) we obtain 

(3-5) J^[QK)] + a M-^[QM] = eQ'(u e )-^[Q(u e )}, o(-) = /'(■)• 

Next we denote 

(3.6) w.(x,t) = ~[Q(u.(x,t))] = Q'(u t )-~[u e (x, ()]. 

Differentiation of (3.5) w.r.t. x yields 


lo n\ dw ‘ , r t \ Q"( u e) , a '( u c ) 2 

(3.7) w + K«.) - *g ?= y-J- S - + = «9 (“.)• jp- 

By a standard regularization argument of the Heaviside function (which we omit), (3.7) 
implies for io+ = H(w e )w e , 

(3.8) M + [„(„,) _ + 0K + )’ < «?(«.)££. 

The maximum principle shows that (3.8) is dominated by Ricatti equation 


— [supu>+(®,*)] + /?[supu;+(x,t)] 2 < 0, 
at x x 


which in turn leads to (3.9), 


(3.9) 


I + i<3(u,(.,i))] = 3Up W'{x, t) < - - — 

X wt (x,0) P 


*+[<?(».(•■ 0))] 

l+jStL+WW-.O)))- 


Finally, we treat the possibly degenerate case where Q' > 0. As in [18], we introduce a 
further regularization where Q{u e ) is replaced by Q s {u e ) = Q(u e ) + 8u e , and (3.4) is then 
recovered by letting 8 J, 0. 


REMARKS. 

1. Theorem 3.1 covers the convective porous medium equation where Q(u e ) = |u e [ 7 it e , 
7 > 0. 


2. In view of the convexity assumption on /, we may set Q(u e ) = a(u e ) and obtain 

(310) ttSSj- 

Since in this case (where Q'(u e ) > a > 0),u e (x,t) converges strongly to the entropy solution, 
u(x,t), of (3.1), the e J, 0 limit of (3.10) also gives us 


(3.11) 


L + [a(u(.,t))] < 


L + [q(u(-,0))] 
l + tL+[a(u(.,0))]’ 
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The a priori estimate (3.11), stated in somewhat weaker form, can be found in [5, Theorem 
1], It shows that the compact solution operator of the nonlinear problem (3.1) tends to 
‘linearize’ the problem as a x (- f t ”f oo) = 0. As noted in [5], the inequality (3.11) requires f 
to be merely C\ and is sharper than both - Oleinik’s result (3.2) as well as its generalization 
in [2, Proposition 1] which apply to C 2 fluxes. 


3. Theorem 3.1 implies, in particular, that the positive variation of Q(u e (-,t )) supported 
on any compact domain is bounded. Consequently, the total variation of Q(u e (',t)) over 
such domains is upper-bounded by 


(3.12) 


||Q(u e (-,f))||j 3 v < Const 0 


£ + [QK(-> Q))] 

l+(3tL+[Q(u e (;0))} 


< 


^ L max <7 L + [a(u e (-, 0))] _ 

< Consto min o' ' 1 + fI+[a(u e (-,0))]’ 


where Consto equals twice the size of supp(?(u e (-,t)) We note that (3.12) holds for BV initial 
data, u e (-,0), over arbitrary domains (with a different coefficient Consto). 


Equipped with Theorem 3.1, we now turn to the main of this section - the local conver- 
gence rate of small viscosity solutions for uniformly parabolic equations 

(3.13) |[«,(«,l)]+^/K(<,l))l = 4 [eW ‘'' ))l 

The difference between u e and its entropy limit u, 


e e (x,t) = u e (x,t) - u(x,t), 


satisfies the error equation 

(3.14a) -Q^[e e (x,t)] + -^[a e (x,t)e e (x,t)] = £-g^[Q(u e (x,t))], t> 0, 

where a e (x } t) denotes the mean- value 

(3.146) a e (x,t)= T a{£u e (x,t) + ( l-O*0MM- 

£ — 0 

We recall that Theorem 3.1 applies to both u e ( x, i) and its entropy limit ii(x,t). In particular, 
in view of the strict convexity of Q, (3.4) implies 


r+r t t> o 
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and in view of the strict convexity of o, (3.11) implies 

£♦[*.«)] 

a 1 -f fL + [a(u(.,0))j 

Inserting the last two inequalities into (3.14b) and using the convexity of a(-,t ) once more, 
we find after little rearrangement the following. 

PROPOSITION 3.2 (OSLO). The averaged convective velocity a e (x,t) given in (8.14b) 
satisfies the OSLC 


(3.15) 


i+ [“. (■,<)] £ ”>(<). m(0 = 


REMARK. The constants Lq < oo and 7) > 1 are given respectively by 
(3.16a) L+ = max(T + [a(u e (.,0))],I + [aK.,0))]) 


(3.166) 


V = 


max a' ■ max Q' 
min a' • min Q' 


Let us now form the dual problem of the error equation (3.14). This is given by the 
backward linear transport equation 

(3.17a) ^[<f> e {x,t)] + a e (x,t)-^-[^ e (x,t)] = 0, t < T, 

with, say C l data, independent of e prescribed at t = T, 

(3-17 b) <f> e ( x ,t = T) = 4>(z). 

Although a e (-, t) may - and in the generic case, will be discontinuous, Proposition 3.2 tells 
us that Theorem 2.2 applies in this case; namely, by (2.27) we have 

(3.18) IIM-.*)l!i*. < 

We are now ready to proceed with our main result announced earlier. 

We convolute the error equation (3.14a) against over its compact support; we 

convolute the adjoint equation (3.17a) against e e (-,t); their sum results in 

j f \ 2 

(3.19) ^e e (-,<) * <f> t (-,t) = e ^2 
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Equation (3.19) governs the evolution of moments of the error. By (3.12) and (3.18), its 
RHS is upper-bounded by 


(3.20) 


eConsto ■ L* ■ 0 < t < T. 

0 max Q (1 +tL£) r >+ 1 "dx u 


Temporal integration of (3.19) together with (3.20) yields 


THEOREM 3.3. Let u(x,t) and u e (x,t ) be the entropy solution and the corresponding 
viscosity solution of (S.l) and (S.1S), respectively. Then there exist constants K > 0 and 
7 } > l 1 such that for any C 1 - function, (f>(x), the follovnng estimate holds 

(3.21) |k(.,TW-<‘(-,TW|| l - < *(1 + W’ +I |I^IU--[£ + ll /%,(•, 0)-«(-,0))|| t .]. 

The dependence of the error estimate (3.21) on the initial data, u(-,0) and u e (-,0), is 
reflected by the following two quantities: 

1. The W -1 * 1 size of the initial error. Of course, in order to obtain the desired 0(e) 
convergence rate, we shall need a rather weak consistency assumption in this direction, 
requiring 

(3.22) ||w e (-,0) — u(-, 0 )||w-m < Const. • e. 


2. The one-sided Lipschitz size of the initial data, measured by 

= max(L + [a(u e (-,0))], L + [a(u(-,0))]). 

We restrict our attention to initial data for which Lq is finite, i.e., 


(3.23) 


ess sup 

x^y 


'u e (x,0) -u e (y,0Y 


x - y 


, ess sup 


Mx,0)-u( y ,0) '| t ^ Const 

V * -y J 


In other words, (3.23) assumes general initial data as long as they do not contain non- 
Lipschitzian increasing discontinuities; in particular, arbitrary C 1 initial data are permitted. 
It would be desirable to extend our result to arbitary BV initial data. In this context, the 
reader is referred to [4] for an almost optimal convergence result for the case of monotonically 
increasing initial data. 

1 The constants K and y may depend on the Btrictly positive extreme values of f" and <f>' (and K may 
depend on the size of supp^), but otherwise are independent of c and t. 
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We have shown 


THEOREM 3.4 ( Uniform estimate on the moments ). Let u(x,t) and u e (x,t) be the en- 
tropy solution and the corresponding viscosity solution of (S.l) and (3. IS) respectively. We 
assume that 

A. The initial viscosity data, u e (-,0), is consistent with the initial entropy data, u(-, 0), 
in the sense that (3.22) holds. 

B. The increasing part of the viscosity and entropy initial data is Lipschitz, i.e., (3.23) 
holds. 

Then for any T > 0, there exists a constant K = K(T) such that for all C 1 -functions, 
4>{x), we have 

(3.24) IW-, TW- u(; T) * ^|U- < K ■ t • II^IU- 

REMARK. Theorem 3.4 suggests W -1,1 as a possible weak topology to study the conver- 
gence of u e . P. D. Lax [7, p. 191] used a similar setup to prove convergence (without a rate 
estimate) already in 1954. 

The estimate on the moments of the error in Theorem 3.4 can be converted into an L 1 
error estimate at the expense of ‘losing’ an additional factor of i/e. This results in the usual 
0(y/e) L 1 -convergence rate in agreement with, [6], [12]. We close this section with 

COROLLARY 3.5. ( L p -error estimate ). Assume that the conditions of Theorem 3-4 hold. 
Then for any T > 0 and p > 1 there exists a constant K (which depends on T and 2=I } but 
otherwise is independent of e) such that 

(3.25) Uu'(;T)-u(;T)lk,<K-e*, p > 1. 

PROOF. Let C(x) be a Cq function with unit mass. For any <j>eL°° we consider 

(3.26) e c (-,T) * <f> = e e (-, T)* fa + e e (-,T) * (<f> — <f>s), <t>s = <t> * 


By Theorem 3.4, the first term on the right does not exceed 

(3.27) ||e,(*,T) * 4>6\\l°° <Ke- ||^||l~ < Const.^ • ||^| |l«. 
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By (3.12), u.(-,T), and likewise, u(-, T ), have bounded variations, and hence the second term 
on the RHS of (3.26) is upper bounded by 

e t (-,T)*(<j>-fa) = [ e e (y,T)fa- -y)dy- 

Jy 

- J v e c {y } T) Jj{--y-z)- s C{- 5 )dzdy = 

= J[ J \e e (y,T) - e t (y - z,T)]fa- - y)dy)±((^)dz < 


(3.28) 


< Const. 5 ||<^||l°° 

Inserting (3.27) and (3.28) into (3.26) we obtain 

(3.29) | J e e (x,T)fax)dx\ < Const.(^ + 6 )||^||l°°. 
Choosing the free parameter 6 ~ y/e, (3.29) with <j> = e^ _1 (-,T) yields 

(3.30) ||e,(.,T)||^ < Const.^|| ee (-,r)||r~ 1 , 
and the result (3.25) follows. 


4. LOCAL ERROR ESTIMATES 

In the previous section we proved that the error between C ,1 -moments of u c and the 
corresponding C 1 -moments of u is of order e. In this section we show that given u e , we can 
use its moments in order to recover pointwise values of u (and any of its derivatives) with 
an error as close to e as desired. 

The idea of such recovery technique is not new. The reader is referred to Mock and Lax 
[9] who post-processed difference approximations in order to recover accurately the pointwise 
values of discontinuous solutions for linear hyperbolic equations. Their approach was later 
extended in [1] and [3] to include spectral and ^dospectral approximations for such equations. 
We shall give a bird’s eye view of both approaches to the post-processing technique in the 
present context of nonlinear equations. 

Assume that u(-,T) is smooth in some fixed ^-neighborhood of x, say 

(4.1) u{;T)eC p {x-e,x + ei p> 1. 

Let fax) be a function supported on the (-1,1) interval such that 

(4.2) J ^fax)dx = 1 , j ^x k fax)- 0, k = 1,2, •••,?— 1 . 
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With 4>s{x ) = |<^(|), where 6e(0, 6) is a free parameter to be chosen later, we have 

(4.3) |(«(-,T) * &)(*) - u(x,T ) | < • IM^)||w*~(*-s,*+fl)- 

The error bound in (4.3) can be made small in two ways: by decreasing 8 - the size of suppifo, 
as was done in [9]; and by increasing p - so that further cancellation due to the oscillatory 
behavior of <f> occurs, as was done in [3]. We shall the sketch the details of both ways. 

By Theorem 3.4 we have 

(4.4) !(«.(-, T) *&)(*) -(«(•, T)*<f> s )(x)\ < K •£■ ll^flk~ 

We now have two types of error estimates: mollification in (4.3) and approximation in (4.4). 
If we choose the free parameter 8 as 

(4.5) ^ = ) ’ l u(p) l loc “ ll u ('» J ’)llw r »-.“>(x-s,x+9), 

then the contribution of the two error terms is of the same order and (4.3), (4.4) yield 

(4.6) |(u e (-,T) * <j>s)(x) — u(x, T)\ < Constp • (1 + lu^lioc)^ • . 

Thus, the ^-moments of u e (-,T) recover the pointwise values of u(-,T) with accuracy as 
close to e as the local smoothness of u(-,T ) permits. 

The estimate (4.6) makes use of a dilated localizer 4>{x) whose support lies in the interval 
(-1,1). If, instead, we choose the support of <f>(x ) to lie in, say (0,1), then (4.6) holds with 
|u^|ioc restricted to [ x+,x + 6), 

(4.6) |u (p) |loc = ||u(-,T')||iyp l0 o[ a . +iI+fl ). 


Thus we axe able to recover pointwise values of u(-,T ) up to the discontinuity. 

In a similar manner, we can recover any spatial derivative of the entropy solution, 
^u(-,T). Indeed, Theorem 3.4 gives us 


fta ft* € 

( 4 ‘ 7 ) - k ' £ ^ h ^^ L00 - K i *+».; 


d*+y 

dx *"*' 1 



Hence, by choosing 


(4.8) 8 = 6 f j-TjT | u (p+»)| — J » N^^lioc — || u (’)^ 1 )l|wp+*'“(i+e»upp<^)) 

we obtain a generalization of (4.6) to higher derivatives, which we summarize in 
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THEOREM 4.1. Assume that the conditions of Theorem 4- 8 hold. Then if we choose <fo(x) 
as described above (4» 2), (4*8), the <f)$-moments of -^u e (' } T) recover the pointwise values of 
% t(-,T) and its derivatives , and the following estimate holds , 

(4.9) |(^u,(-,T) * *)(*) - ^«(*.r)l s Constp • (1 + • E5* 1 , 

Next, let us consider the generic case of smooth initial data, tt(x, 0)e5. Then, the entropy 
solution of (3.1) (with C°° convex flux /) is piecewise smooth [13]. In this case, we are able 
to recover the pointwise values of -^u(x, T ) with error as close to e as desired, if we take p 
large enough in (4.9). 

We close this section with a brief description of the spectral post-processing technique 
[3] which enables the pointwise recovery of u(-,T) and its derivatives. 

We restrict our attention to the case (4.1) where the symmetric interval (x + 6, x + 6) is 
free of discontinuities of u(-, T). Let £(x) he a Co°( — 1> 1) function, normalized such that 

(4.10) C(z = 0) = 1, 


and we set 


(4.11) 


i i t 1 'i x \n I x \ n /' ) - 1 sin ( w + 5)£ 

M*) = D «((> = Jin I f 


We note that in this case the support of the regularization kernel, <^/-(x) is kept fixed. 
Instead, by increasing N , we obtain a highly oscillatory kernel whose monomial moments 
satisfy (4.2) modulo a spectrally small negligible error. 

Standard error estimates for the truncated Fourier projection, Sjf, give us [3, Section 3] 

IK-, T) * - u(x, T ) I = |(/ - S N )[u(x - 6; T)C(-) - u(x,T)m\ 

( 412 ) i n ft 

< Constp • -^|K-,T)||H^.°°(*-«,*+fi)- 

Applying Theorem (3.4) with <f> = <f>N we find 

(4.13) |(u,(-,T) * 4>n)(x) - ( u(-,T ) * ^Ar)(a=)| < K e - ^ Const.e • N 2 , 


Hence, by choosing 
(4.14) 1 




,«l, - 


loc = ||u(-,r)||lVf-“(*-e,x+fi) 


we recover from (4.12), (4.13) the local error estimate we had in (4.6), namely we have 


(4.15) 


l(u.(-,rWw)(»)-«<x.T)| < Const, >WC ■£*■ 
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5. CONCLUDING REMARKS 


The results of Section 4 hinge on the a priori estimate of the moments in Theorem 
3.4, which, in turn, is based on the W 1,0o -energy estimates for the OSLO satisfying linear 
transport equations studied in Section 2. In this section we provide still another derivation 
of Theorem 3.4 which amplify the direct linkage between the OSLO and Theorem 3.4. To 
demonstrate our point we will concentrate on the, say 27T, periodic case. 

Let 


(5.1) 


E e (x,t) = J u e ((,t)d( - J u{(,t)d( 


denotes the difference between the primitives of the small viscosity approximation, U e (x,t ) = 
/*«,(£, t)d£, and that of the entropy solution, U(x,t ) = / x u(£,t)d£. Integration of the error 
equation (3.14a) gives us 


(5.2) 


— [E e (x,t)] + a e (x,t)-Q^[E e (x,t)] = e-?j-[Q(u e (x, i))]. 


We shall now energy estimate (5.2) along the lines of our study of irreversible linear transport 
equations in Section 2, compare (2.16). Integrating (5.2) against sgnE ei we obtain 


(5.3) 


£ 

dt 


f d 

\\Ec{-,i)\\v + Jd e (x,t)—[E e (x,t)]dx = 


f d 

= £ J sgnE e (x,t)— [Q(u e (x } t)))dx. 


Integration by parts of the second term on the RHS of (5.3), together with the upper bound 
of the RHS by the BV estimate of Q(u e (-,t)) in (3.12), yield 


(5.4) 


d f d 

—\\E e (-,t)\\ L i < J — [a e (x,t)} ■ E e (x,t)dx + Const. e. 


Thanks to the judicious positive sign on the right, we may now use the OSLO (3.15) to 
obtain 

d 


(5.5) 


dt 


!!■£«(•, *)IU l ^ m(t)||E«(-,<)|| t i + Const .e. 


By integration of the last inequality we conclude - in agreeement with Theorem 3.4, that 


(5.6) 


U.(;t)-U(;t)\\ L i <«•■£, 0<f<r. 


Using this we can derive L°° (and consequently, L p ) estimate as follows. Let E e (x, t) assumes 
the Fourier expansion 

(5.7) E.(x,t)= ■£ E(k,t)e “ «. 

|fc|<oo 
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By (5.6) we have 

(5.8) |£» e (fc,t) < 2^||-®«( , iOIU l — 2jt e ‘ 

Moreover, since both u e (-,t) = '§^.U e {- ) t') and n(-,t) = are uniformly BV (w.r.t. e), 

their Jfc-th Fourier coefficients decay like ^ |fc| -1 , and hence 

(5.9) \E e {k,t)\ < Const. -|p |Jb| > 1. 

Inserting (5.8), (5.9) into (5.7) we find 

(5.10) ||£,(',<)|!l- < E £•« + E Coost.rij < Const. 

1*15* W>* 1 1 

Standard interpolation between the If error estimate in (5.6) and the L°° estimate in (5.10) 
gives us 

THEOREM 5.1. Assume that the conditions of Theorem 3-4 hold. Then for any T > 0 
there exists a constant K = K(T) such that 

(5.11) \\U t (;t)-U(;t)\\^>K-e^r, p > 1, 0 < t < T. 
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